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We study constraint effective potentials for various strongly interacting 0^ the- 
ories. Renormalization group (RG) equations for these quantities are discussed and 
a heuristic development of a commonly used RG approximation is presented which 
stresses the relationships among the loop expansion, the Schwinger-Dyson method 
and the renormalization group approach. We extend the standard RG treatment 
to account explicitly for finite lattice effects. Constraint effective potentials are then 
evaluated using Monte Carlo (MC) techniques and careful comparisons are made with 
RG calculations. Explicit treatment of finite lattice effects is found to be essential in 
achieving quantitative agreement with the MC effective potentials. Excellent agree- 
ment is demonstrated for d = 3 and d = 4, 0(1) and 0(2) cases in both symmetric 
and broken phases. 
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I. INTRODUCTION 



Strong-coupled relativistic quantum field theories (RQFT's) play an essential role in any 
modern understanding of subatomic phenomena from the Higgs mechanism to fundamental 
theories of nucleon and nuclear phenomena. There are important theoretical issues, such 
as spontaneous symmetry breaking, which can be studied with relatively simple scalar field 
theories. Unfortunately, even the simplest interacting relativistic quantum field theories are 
notoriously difficult to solve. A classic example is the issue of triviality in 4-dimensional 0^ 
scalar field theory which has been tackled over the years using nearly the full complement of 
weapons in the theorists' arsenal (loop expansion, strong coupling- and e-expansions, trun- 
cated systems of Schwinger-Dyson (SD) equations, renormalization group (RG) techniques 
and Monte Carlo studies), without yielding a definitive answer. (See, e.g., Ref. .) It is 
self-evident that accurate and relatively simple approximate methods for solving RQFT's 
would be of enormous value. 

In the present paper we discuss approximations to (f)'^ theories emphasizing throughout 
comparison with lattice Monte Carlo calculations which, in a sense to be made explicit below, 
we consider exact. In particular we will focus on the uniform properties which can be ex- 
pressed in terms of "effective potentials" (or more precisely, "constraint effective potentials" 
as defined, e.g., in Ref. ) from which, for example, zero-momentum n-point couplings 
may be extracted. The present results may be summarized most succinctly by stating that 
the renormalization group (RG) effective potential - evaluated so as to account for lattice 
artifacts - provides a spectacularly accurate description of its lattice counterpart with little 
calculational complexity for a wide variety of scalar field theories including symmetric and 
broken phases in 2, 3, and 4-D with 0(1) and 0(2) symmetries. In contrast, other meth- 
ods such as the use of truncated SD equations yield only qualitative agreement at best. 
Moreover, as will be emphasized below, the simple RG effective potential used here is known 
to be the first term in a momentum expansion meaning that systematic improvements in the 
approximation are possible. Perhaps even more important is the evidence that - in contrast 
to, e.g., the hierarchy of SD equations - this expansion appears to be strongly convergent 
for cases with strong coupling. Indeed, estimates of next-to-leading order correction to be 
made below are quite small but consistently improve the already extraordinary agreement 
with lattice results. Finally, as will also be addressed below, the RG approach is ideally 
suited to the construction of effective low energy theories which can economically account 
for the long-wavelength behavior of the full theories. We will also discuss some anticipated 
exploitations of this possibility. 

In section II we review briefly the history of the version of the RG effective potential 
used in the present work providing highlights and references to several other recent studies. 
In section III we provide a heuristic derivation of the RG equation for the effective potential 
and compare with truncated Schwinger-Dyson. In section IV we describe the methods we 
use in extracting the effective potential from Monte Carlo simulations and in solving the 
RG equations numerically. In section V we present the results of our calculations comparing 
Monte Carlo and RG descriptions for broken and symmetric phases of 0(1) and 0(2) 0^ 
field theories in 3- and 4-dimensions. 
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II. REVIEW OF THE RG EFFECTIVE POTENTIAL 



We now trace the development of the form of the RG equations used in the present 
study. The original work by Wilson on the renormalization group was in essence a 
discussion of how to obtain the effective action at long wavelengths from a bare action at 
short wavelengths. He and his coworkers discussed a variety of approximations to the full 
(and intractable) RG equations. Shortly thereafter, Wegner and Houghton derived an 
exact RG equation for the effective potential by integrating over thin shells in momentum 
space, a technique referred to as the "sharp cutoff method." (See Ref. p for an early 
discussion of this and alternative "smooth cutoff" methods. Also see Ref. for early 
discussions of many RG concepts. Ref. @] provides a later discussion of some of these 
ideas in the course of constructing a novel and intuitively appealing RG-based proof of 
the renormalizability of 0^ field theories.) A seminal early application of the Wegner and 
Houghton method to 0^ theories was presented by Hasenfratz and Hasenfratz P . (See 



also Ref. |]T0| for an illuminating comparison of analyses of the cutoff dependence of the 



Higgs mass using the Wilson-based and Wegner /Houghton-based RG methods.) These 
authors obtained a tractable approximation to the exact Wegner and Houghton equation by 
projecting the full effective action onto an effective potential for which all n-point couplings 
are momentum independent. They used the resulting RG equation to extract information 
about RG flows in the vicinity of Wilson and Gaussian fixed points in c? = 3 and to establish 
the absence of a Wilson fixed point in d = 4. They also extracted various critical exponents 
for the d = 3 case and found them to be in good agreement with values determined by a 
variety of other techniques. 

There has been a recent resurgence of interest in RG effective potentials. Wetterich 



and coworkers |]1T| , for example, have extended the approach of Hasenfratz and Hasenfratz 
using the so-called "smooth cutoff" method in deriving their RG equations. Their technique 
permits evaluation of, e.g., the anomalous dimension, t], which is implicitly zero in the 
Hasenfratz and Hasenfratz approximation and yields values of critical exponents which are 
in impressive agreement with those determined by other methods. (See Ref. |12| for an early 
calculation of t] using an extension of one of Wilson's original approximate RG equations.) 
This increase in calculational power comes, however, at the cost of a substantial increase in 



numerical complexity. In other work, Morris |TB| gives an in-depth discussion of the RG 
effective potential as derived using both sharp and smooth cutoffs. He also developes a 
version of the systematic momentum expansion alluded to above and presents an analysis 
of 0^ in d = 3 using the first correction to what is in effect the Hasenfratz and Hasenfratz 
method (which appears as the zero-momentum limit in both the sharp- and smooth-cutoff 
derivations). Without attempting to be comprehensive, we mention other related papers 



which appear here as Refs. |15] , [It] and 



III. A HEURISTIC DEVELOPMENT OF THE RG EQUATION FOR THE 

EFFECTIVE POTENTIAL 

Derivations of RG equations for the effective action or their reduction to expressions 
governing the effective potential appear in Refs. P,[7|,P,p|, pT| , pr3| , p!7[] and will not be reproduced 
here. Instead we provide a heuristic and, we hope, intuitive development of the RG equation 
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in the zero-momentum limit which has as its point of departure the famihar expression for 
the one-loop contribution to the effective action of a uniform system. 

In d Euclidean dimensions, the bare action for a one-component 0^ theory is given by 



where 



(1) 



f/o(0) = ^ml<p' + 1ao0^ (2) 

The one-loop perturbative correction to the bare potential for a uniform system corresponds 
to the familiar expression (primes denote differentiation with respect to 0) 

UM = lJ0^,Hk' + U'^m (3) 

which can be interpreted physically as the sum of zero-point energies of vacuum modes for a 
free scalar meson theory with meson mass-squared = Uq (0). This quantity is divergent 
and we regularize it by adopting a UV cutoff Aq: 

UM = lJ-0^, ^(Ao - 1^1) Hk' + U'om (4) 

We note that the lowest order truncated Schwinger-Dyson equation can be obtained from 
the following expression for what we will call the Schwinger-Dyson effective potential, Usd, 
in which the 1-loop contribution t/i(0) as given above is self-consistent in the following sense: 

Usd{<I>) ^ f/o(0) + lJ-0yi ^(Ao - 1^1) Hk' + U'snm. (5) 

The RG equation for the effective potential Urg may be obtained heuristically by starting 
with the loop expansion of the effective action, the first term of which is given in Eq. ^. 
Assume we know f/(0. A), the effective potential at momentum scale A. We then have the 
following expression for U{(f),A — dA), the effective potential at a slightly smaller momentum 
scale: 

f/(0,A-rfA) = f/(0,A) + ^|^0(A-|fc|) e{\k\-A + dA) \n[e + U'^{<P)] 



{2ny 

U{(j), A) + ^A''-^ In [A^ + f/"(0. A)] dA + 0{dA^) (6) 



where = J dVld/{2T[Y = l/{2^~^Tc''-^'^T{d/2) ) and the higher order loops contribute at 
0((iA^). Taylor expanding the left hand side and taking the limit rfA — gives the RG 
equation: 

;?^^?£iM)=_^A^-M„|A^ + f/;„(^.A)]. (7) 

The essential point is that where the Eq. ^ requires Usd be made self-consistent across all 
momentum scales at once, Eq. ^ ensures that Urg is self-consistent at every scale. Put 
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another way, the evolution (or "flow" ) of Urc at some scale A is determined self-consistently 
by Urc at that same scale. We make think of Usd as the result of the crudest possible 
discretization of the flow which yields Urg, namely one in which the integration of the 
differential equation in Eq. |^ is carried out in a single step. 

The extension to 0{N) is straight-forward. Deflning p = \[$^ and with primes denoting 
differentiation with respect to p, we flnd 



dURG{p,^) A 



dA 



-f A'"' {ln[A^ + f/L(P, A)] + (AT - 1) ln[A2 + U'^^iP. A)/p]}. (8) 



To make contact with, e.g., Hasenfratz and Hasenfratz |^, we scale all dimensionful quan- 
tities by the appropriate power of A so as to render them dimensionless. For example, 
— > X = ^/A"'/^^^ (where in this approximation the anomalous dimension, 77, vanishes). 
Designating the rescaled potential as 

URG{x,t)=URG{<P,^)l^^ (9) 

where t = ln(Ao/A), the one-component equation (Eq. ^ becomes, after dropping irrelevant 
constants, becomes 

_ ^^ dURG{x,A-) _ dURG{x,t) _ 
dA ~ dt ~ 

^ ln[l + U'i^ai^, t)] + (1 - d/2) X [/^^(x, t) + d Urg{x, t) (10) 

where primes now denote differentiation with respect to x. This is Eq. 22 of Ref. (with 
rj = 0). Finally, deflning f{x,t) = Urg^x^I), Eq. ^ becomes 

A t" 

f = YiT7 + + (1 + ^/2)/ (11) 

where the "dot" means differentiation with respect to t. This is Eq. 9 of Ref. ^ and is a 
form of the RG equation for the effective potential frequently encountered in the literature. 

These RG equations may be integrated to flnd the effective potential at momentum scale 
A subject to the boundary condition that the bare and effective potentials coincide at the 
UV cutoff scale Aq. In so doing, quantum corrections associated with single meson modes 
with momenta between A and Aq are (approximately) included in the effective potential 
while contributions from modes at lower frequencies must still be treated explicitly. As 
touched upon above, the power of the RG approach derives from its self-consistency at each 
momentum scale. This feature is evident in Eq. |^ where it is clear that the evolution of 
Urg at momentum scale A depends on a quantum correction which is consistent with the 
potential itself at that particular momentum scale. 

Another important and perhaps surprising aspect of the RG equations is that - while 
they have the form of a truncated loop expansion - no such truncation has been made. As 
originally shown by Wegner and Houghton higher-order loops are suppressed by powers 
of the differential dA and their equation therefore has the form of a one-loop result but is 
in fact exact. Thus the RG equation for the effective potential, obtained by Hasenfratz and 
Hasenfratz as an approximation to the Wegner and Houghton expression has nothing to do 
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with explicitly dropping high order quantum corrections but instead depends on projecting 
the full effective action onto a form with independent n-particle couplings. (This is what 
is meant by the "zero- momentum" RG effective potential.) Of course, just as is the case 
when solving truncated SD equations , quantum corrections of arbitrarily high order in a 
perturbative loop expansion are automatically generated in the course of integrating the RG 
equation. (The relationship between the RG equation and the loop expansion was touched 
upon by Polchinski ||^ and addressed in more detail by Morris []TB[ who also made a direct 
comparisons with the SD expression given below in Eq. |15[) 

Once the effective potential is obtained, the renormalized mass and 4-point coupling are 
given by: 



(fa 



m 



(12) 



A 



1 d'Ueff 

6 dd)^ 



(13) 



where 0o is the value of the field at the minimum of f/e//- For example in the case of the 
Schwinger-Dyson equation, Eq. H, in the symmetric phase where 4>o = 0, we have 



m 



SD 



ml + 3\sD J 



(14) 



SD 



^SD and Usoi 



where we have used Ugj^{(f) = 0) 
by Ao yields the familiar (see, e.g., Eq. 2.17 of Ref. 
"dressed" mass: 



= 0) ^ 6XsD- Approximating Xsd 
SD self-consistency relation for the 



m 



SD 



3An 



d'^k 



1 



{2^Y k^ + m 



(15) 



SD 



where the self-consistent nature of the SD method is apparent. (We note that the approxi- 
mation defined by this equation for msD also known [0] as the "cactus" approximation.) 
Likewise, within the same approximation, the 4-point is 



A 



SD 



Ao/ 



l + 9Ao J 



d'^k 



{27iY {k^ + mlj^f 



(16) 



Examination of RG renormalized masses and 4-point couplings following from Eq. ^ (or 
the equivalent expressions given subsequently) requires finding numerical solutions to the 
non-linear partial differential RG equations. In the next section we discuss our numerical 
methods for solving the RG equations and for determining the "exact" effective potentials 
from our Monte Carlo calculations. 



IV. NUMERICAL METHODS 

We have carried out RG calculations for the effective constraint potential in two com- 
plementary ways. In the first, appropriate to studying the continuum limit, we treat A as 
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a continuous variable. In the second approach, appropriate to direct comparison with our 
Monte Carlo calculations which of course refer to a finite lattice, we treat the lattice normal 
modes explicitly. We now give detailed descriptions of both methods. 

In the first approach, we solve Eq. ^by treating A as a continuous variable. (We have also 
evaluated the equivalent expressions, Eq. |T^ or Eq. 0. Although one version is often better 
suited to a given calculation than the others, we found them to be entirely equivalent to one 
another as expected.) We begin by expanding Urg{4>, ^) in a power series in (j) truncated at 

M 1 

[/Ha(0,A) = E 7^l^2^i^^''. (17) 

1=1 

The boundary condition referred to in Section III is imposed by setting 

W2(Ao) = ml W4(Ao) = Ao, and W2^(Ao) = for 2i > 6. (18) 
The range of the A integration is divided into intervals of equal length AA: 

A„ = Ao - nAA (19) 

with 

AA = (Ao-A,^)/Ar, (20) 

where A/^ = 2ti /Niat is the infrared scale appropriate to a lattice of Nf^^ sites and Aq is 
chosen to yield the same phase space volume as the finite sum-over-modes. 

A similar decomposition is made of a range of 0, the other independent variable: 

00 < • • • 0i < 0i+i < . . . < 0j. (21) 

Each step in the integration is accomplished by first making a least squares fit of the poly- 
nomial on the r.h.s of Eq. |l^ to UjiG{(p,-^n) in order to determine the W2i(A„) and then 
computing 

Unai^j, An+i) = Unci^j, K) + AA ■ A^^ ^ \n[Al + U'i^aiK A„)] (22) 

for each j where U'^Q{(f)j, A„) is computed analytically from the polynomial fit. Starting with 
n = 0, this procedure is iterated until ultimately the J +1 values of URG{(pj, ^m) are found. 
We refer to this as the "continuum RG effective potential" (C-RG). Obviously this method 
yields directly the m-point couplings, Wm(A/i?,), as well. The calculation can sometimes be 
unstable in the vicinity of the 0-endpoint, 0j, and some care is required in choosing values 
of M, N and J to ensure stable integration. In the calculations to be presented below, we 
typically used 

M ~ 10, ~ 4000, and J ~ 50. (23) 

In our second approach, an alternate RG calculation was formulated to facilitate com- 
parison with the lattice results to be described below. The derivative terms in the action 
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(Eq. H) are approximated on the lattice by a finite difference formula. For a uniform sys- 
tem on a Nf^^ lattice with periodic boundary conditions, the vacuum modes are momentum 
eigenstates with wavenumbers in each direction given by k„ = 2mi/Niat, n = 1,2, . . . Niat 
where we take the lattice constant to be unity, a = 1. We make the following identification: 

4^.f(^')-Wr ^ ^•••^ fikl + kl + ... + kl) (24) 

l^^J ^^lat ni=l n2=l n^=l 

where fc^ = 4sin^(K„/2). Now define the A*";"^^ scalar quantities 

then sort them by magnitude and relabel: 

kl>...kl>kl^^>...>kl, (26) 

lat 

where kl = M and k'^d = 0. Our "latticized" version of the RG equation for the effective 

lat 

potential, Eq. is 

URG{<f>^, fc^+i) = UnG{<f>^. k'^) + ^ ^ ln[fc^ + U';,ai<l>^. klf)] (27) 
with the boundary condition 

f/«G(0.,A;D = f/o(0^). (28) 

We will refer to this as the "lattice RG effective potential" (L-RG). 

The lattice Monte Carlo calculations to be presented below employed (i-dimensional cubic 
lattices with periodic boundary conditions in conjunction with the Metropolis algorithm to 
sample the 0-configurations with weights given by the exponentiated action. The Monte 
Carlo observable of primary interest here is the "constraint effective potential" which we 
refer to as the "MC effective potential" and designate by t/e//- Formally, this quantity is 
defined via 

exp[-<f/,^;(0)] = I V{<^} 5(0 - exp[-^[0]] (29) 

i 

where X^i means a sum over all lattice sites and is the average field on the lattice for a 
given field configuration in the path integral. 

Since the Monte Carlo method generates sequences of field configurations, weighted 
by exp[— ^[^jn], f/e//(<^) is extracted from the lattice Monte Carlo calculation from a his- 
togram of the number of configurations generated vs the average value of the field for those 
configurations. Let dMi^^ be the number of configurations with average field values in an 
interval c?<^ about 0. Then 

dM{^) oc exp[-iV,l,f/e//(0)] rf0 (30) 
and, to within an overall constant. 
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In 



(31) 



For 0{N) theories this expression is modified shghtly. Let p be the average value of p = 

on the lattice and let dj\f{p) be the number of configurations in an interval dp about 
p. Then, since 



dM{p) (X p^-' exp[-iVf,, U^ffip)] dp, 



we have 



Ueffip) 



In 



1 dU{p) 



p^-i dp 



(32) 



(33) 



Of course the most probable average field values are near the minimum of f^e//- In order 
to determine f/e// away from its minimum; i.e., to sample a range of relatively improbable 
values of (or p), it is often necessary to add a constraint to the updating algorithm which 

m 



rejects "moves" which take the average field value outside a preset "window'' 



By 



matching Ueff^s for a number of adjacent windows, it is possible to detemine the lattice 



effective potential over a wide range of average field values. (See Ref. [jT9| for an interesting 
alternative method.) 

Finally we comment on the "convexity" of U^ff. It can be shown by very general argu- 
ments that the condition U^j-j^cp) > (primes denote differentiation with respect to 0) 
must be satisfied for all (p. Presently, we will show Uef/s for "broken" symmetry phases of (p^ 
theories which do not satisfy the convexity condition. The apparent contradiction is resolved 
by noting that convexity is required only in the thermodynamic limit; i.e., as Ni^t —>■ oo. 
For the finite lattice results to be presented below, convexity is not a fundamental property 

of Ueff. 



V. COMPARISON OF MC AND RG EFFECTIVE POTENTIALS 

We now present MC effective potentials as defined in Eqs. ^ or |53| for a variety of cases. 
We also show the corresponding continuum and lattice RG effective potentials of Eqs. ^ and 

respectively. In computing the former RG potential, we always set the IR cutoff; i.e., A^r 
in Eq. |T9|, to 2'7i/Niat- Fixing the UV cutoff, Aq, is ambiguous. We adopt the prescription 
that Aq should be chosen so that the phase space volumes integrated over should be the 
same for continuum and lattice calculations. This correspondance is made quantitative by 
setting /(A;^) = 1 in Eq. |2^ which implies, in effect, 

A, / k''-Uk = A,^ = l. (34) 
JO d 

Then, e.g.. 



Ao{d 


= 2) 


= (47r)V2 


Ao{d 


= 3) 


= (67r2)i/3 


Ao(rf 


= 4) 


= (327r^)i/^ 
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In the calculations shown in this work the bare coupling was fixed at Aq = 10 while the bare 
mass mo was adjusted to give a desired renormalized mass. The MC effective potentials 
were fit with a polynomial of the form 



Ueff{<P) = lm%<p' + ^Xfu<P' (36) 

to determine the renormalized MC masses and 4-point couplings. 

In our RG calculations we also kept the 4-point coupling fixed at Aq = 10. We now 
describe our procedure for determining the bare mass in the RG calculations. For symmetric 
phases, the RG bare masses were adjusted to yield d'^URc/d(p'^\(j,=o = rn'j-^. For broken 
phases, the RG bare masses were adjusted so as to give the best fit "by eye" to Uejf in the 
vicinity of its minimum (see discussion below). Since the MC and RG effective potentials 
in the symmetric phase behave identically at small by construction, we typically compare 
- for both symmetric and broken phases - MC and RG results for the subtracted quantity 
U — ^m'j-^cj)'^ rather than results for the ?7's themselves. Table 1 summarizes results of 0(1) 
and 0(2) calculations in both symmetric (S) and broken (B) phases in d = 3 and d = 4. For 
the RG results, m = a/Z^ and A = W4/6. The Table shows that our bare masses squared are 
typically ~ —2^ while the renormalized masses squared are ~ ±0.2^ to ±0.3^ with 
renormalized couplings of A ~ 1 to 4. 

Clearly quantum corrections are large and simple perturbation theory {e.g., the loop 
expansion) cannot account for them. For comparison purposes we examine results from the 
simplest Schwinger-Dyson (SD) equations ||^. Consider, for example, a "latticized" d = 3, 
0(1) version of Eq. [15| : 

1 lat 1 



= ml + 3Xo—^ E p , ^2 (37) 

^^lat J\f=l f^Af + ^SD 

where the fc^'s are defined in Eq. In this expression, a value of = —2.61^ is required 
to give msD = 0.285 for Aq = 10 and Niat = 16 while the corresponding MC calculation 
yielding mfn = 0.285 requires = —2.16^. The latticized expression for the SD 4-point 
coupling, Eq. ITS, is 



A5D — Aq/ 



lat 



1 + 9^0 ^ E , ™2 )2 ' (38) 

^^lat J\f=l y'^N + ^SD) 



and gives A^^i = 0.65 as compared to the MC value of 1.49. As another example, we 
find that, for d = 4, 0(1), the MC-determined renormalized mass is mfn = 0.276 using 
ml = —1.88^ compared to the latticized SD value ml = —2.13^ needed to yield the same 
renormalized mass. Similarly, the latticized SD equation for the 4-point coupling gives 
XsD = 1-28 vs Xmc = 3.42. Not surprisingly, these simple SD equations do not in general 
provide a quantitatively accurate picture of the quantum effects. Although we do not display 
SD results for the effective potential, it is apparent from the 4-point couplings that the SD- 
derived effective potential will yield at best only a qualitative description of the MC effective 
potential. 

We turn now to the Monte Carlo calculations of the effective potential and comparisons 
with the two RG effective potentials, C-RG and L-RG, discussed in the previous section. 
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As mentioned previously, since we adjust the bare mass in the RG calculations to give the 
MC fit mass, mfn, in most of the figures that follow we display the subtracted quantity, 
U — ^m'j^^(j)'^, to emphasize the predictive power of our calculations. 

Figure |l| displays the subtracted {i.e., U — |mj^^0^) MC effective potential for d = 3, 
0(1) in the symmetric phase on a 16'^ lattice. A MC bare mass squared of —2.16^ yields 
^fit = 0.285 and Xfu = 1.49, as mentioned above and shown in Table 1. Note that we 
have used the "windowing" technique described in Section IV to extend the range of at 
which the MC effective potential is determined from a maximum value of about 0.1 which 
is feasible without windowing to a maximum value of 0.22. (The relative normalizations of 
adjacent windows were determined by matching intermediate values well determined in both 
windows.) However, the values of mfu and Xfu were determined by fits to the MC values of 
Ueff for < 0.1 only. To reproduce this value of mfu, a C-RG bare mass squared of —1.91^ 
was required. Of course this numerical value depends on our prescription for Aq which was 
discussed above. (However, the C-RG effective potential is only weakly dependent on Aq 
provided the RG bare mass is adjusted to give a particular value of the renormalized mass, 
m.) The C-RG value for the effective 4-point coupling is A = 1.89 and it is evident from 
Figure that the C-RG calculation (solid line) does not accurately reproduce the quartic 
and higher 0-dependence of the MC effective potential. By comparison, a bare mass squared 
of —2.14^ was required in the L-RG calculation to give m = mfu. That the L-RG bare mass 
differs from the MC value by less than 1% is especially gratifying; in contrast, the C-RG 
bare mass differs by about 12% while the SD difference is also about 12% but in the opposite 
direction. The L-RG value of A = 1.43 is also quite close to Xfu = 1.49. However, inspection 
of Figure |I| — where the agreement between MC and L-RG results over the wide range of 
spanned using the windowing technique is spectacular — suggests that even this small 
discrepancy is probably due to the fact that the lattice Urq contains terms of order 0^ and 
higher while our fit to Ueff does not. We note that our limited MC "data" are not good 
enough to extract the 0^ or higher terms reliably. Also appearing in Figure |^ as the dashed 
line is the L-RG result omitting all contributions of order 0^ and higher. Comparison of 
this curve with the MC results show that these high order contributions are important for 
large values of and that the L-RG accounts for them quite well. To summarize, the L-RG 
calculation which includes lattice effects in the integration of the differential equation for 
Urg yields excellent quantitative agreement with MC effective potentials even on a relatively 
small lattice provided a very small adjustment in the RG bare mass is permitted. For the 
C-RG calculation, comparison with the MC results is not straightforward because of the 
ambiguity in choosing Aq which then results in ambigiuties in the bare mass. Futhermore, 
even when the bare mass is allowed to vary freely, agreement is poorer than for the L-RG. 

Unsubtracted effective potentials for the d = 3, 0(1) case in the broken phase appear in 
Figure ||. As mentioned above, the bare masses were adjusted in the RG calculations to yield 
the correct location of the minimum. From Figure ^ we see again that the description of 
the effective potential provided by the L-RG calculation is excellent. The C-RG calculation, 
however, is incapable of simultaneously reproducing the location in and the depth of the 
minimum in the MC effective potential. The fit shown is a (non-unique) compromise. The 
subtracted potentials for this case are shown in Figure |^ and again demonstrate the superior 
description of 0^ and higher contributions to the effective potential provided by the L-RG. 
Indeed, the description of the MC results by the L-RG is again spectacular. The bare mass 
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squared used in this calculation is —2.23^ which differs from the MC value of —2.24^ by less 
than 1%. (The high quality of this fit again leads us to attribute the differences between 
Xfit and A from the L-RG calculation — 1.48 and 0.52, respectively — to the absence of 0^ 
and higher terms in the fitting function rather than to shortcomings of the L-RG method.) 

Similar findings were obtained for d = 4, 0(1) as well as 0(2) cases in ci = 3 and 
d = 4 in both symmetric and broken phases. Values of the bare and renormalized masses 
and 4-point couplings are reported in Table I. The L-RG description of the subtracted MC 
effective potential is nearly perfect in all cases and is distinctly superior to that provided 
by the C-RG. We note that this is true even when the correlation length {i.e., m~^) is large 
compared to the overall lattice size. For our d = 4, 0(2) case, a small MC mass of 0.06 was 
found {i.e., the correlation length was ~ 15) on an 8^ lattice (!) and, yet, as shown in Figure 
^, the L-RG calculation agrees almost prefectly with the subtracted MC effective potential. 
Furthermore, the L-RG bare masses never differ by more than 1.5% from the MC values 
and are typically much closer. Overall, we conclude that the L-RG method provides an 
extraordinarily accurate, semianalytic description of lattice MC calculations of the effective 
potential at a tiny fraction of the computational cost. 



VI. SUMMARY AND CONCLUSIONS 

We have presented a discussion of the renormalization group (RG) effective potential 
for 0^ scalar field theories which stresses its relationship to the loop expansion and the 
Schwinger-Dyson (SD) method. We show in a simple way that the unique feature of the RG 
method - which is presumably the origin of its ability to describe the full theory as well as 
it does - is its self-consistency at all momentum scales. We have formulated a version of the 
RG equation for the effective potential which explicitly treats finite lattice effects so that RG 
calculations may be compared directly with lattice Monte Carlo (MC) calculations without 



extrapolations of the latter to infinite lattice size and zero lattice spacing (See Ref. and 
references contained therein for some examples of such extrapolations). We compare such 
"latticized" or L-RG potentials along with standard continuum or C-RG potentials to 
constraint effective potentials obtained via MC methods. The L-RG effective potentials are 
in excellent quantitative agreement with the MC results for d = 3 and d = 4, 0(1) and 
0(2) cases in both symmetric and broken phases. More specifically, we find that, upon 
adjusting the L-RG bare masses so as to yield renormalized masses equal to those found 
in the MC calculations, the L-RG bare masses never differ by more than 1.5% from the 
MC values and are typically much closer. The renormalized 4-point couplings found by 
solving the L-RG equations are also typically in close agreement with the MC values. In at 
least one case for which the MC effective potential was determined over a wide range of (j) 
using windowing techniques, there is strong evidence that 0^ and higher contributions are 
also well accounted for in the L-RG approach. In another case, we showed that the L-RG 
approach works very well even when the correlation length is appreciably larger than the 
overall lattice size (m Niat — 0.5). Comparison of C-RG calculations with the MC results 
was found to be ambiguous; however, even when ambiguous parameters were tuned to give 
optimal agreement with the MC quantities, this accord was found to be inferior to that 
achieved using the L-RG without tuning. Finally, latticized SD equations were seen to yield 
only a rough qualitative description of of MC masses and 4-point couplings. 
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We conclude that the L-RG method provides an very accurate, relatively simple descrip- 
tion of lattice MC calculations. It may prove to be a useful guide in extrapolating finite 
lattice results to the continuum limit since the relation between the L-RG and the C-RG 
is simple and well-defined. We may also be able to apply it to some specific problems of 
current interest such as accounting for the "Universal Effective Potential" for 0^ field the- 
ories in (i = 3 recently obtained via MC techniques by Tsypin |jT9|. The success of the 
method encourages us to make various extensions of it. For example, we plan to go beyond 
the "zero-momentum" approximation so as to be able to compute wavefunction renormal- 
izations ]T^pil|] . We have also begun to examine including Yukawa-coupled fermions in the 
approach which would be a crucial step toward our long term goal of constructing effective 
theories of low-energy nuclear phenomena whose field theoretic structure is fully understood 
and for which we have an accurate and relatively simple approximation scheme of proven 
reliability. 
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FIGURES 

FIG. 1. Subtracted effective potentials for the d = 3, 0(1) case in the symmetric phase. The 
solid dots are the MC results while the solid and dotted lines correspond to the C-RG and L-RG 
calculations, respectively. Statistical uncertainties for the MC effective potential are not shown 
but are typically smaller than the size of the dots. The dashed curve is the L-RG result with cp^ 
and higher contributions omitted. Input and output quantities are discussed in the text and are 
presented in Table I. 

FIG. 2. Effective potentials for the d = 3, 0(1) case in the broken phase. The solid dots are 
the MC results while the solid and dotted lines correspond to the C-RG and L-RG calculations, 
respectively. Input and output quantities are discussed in the text and are presented in Table I. 

FIG. 3. Same as Figure |2| but for the subtracted effective potentials. 

FIG. 4. Same as for Figure |3| but for d = 4, 0(2) in the symmetric phase. 
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TABLES 



TABLE I. Input and ouput quantities for Monte Carlo (MC) as well as continuum (C-RG) and 
lattice (L-RG) renormalization group calculations for d = 3 and cZ = 4, 0(1) and 0(2) cases in 
both symmetric (S) and broken (B) phases. The bare masses used are denoted by mo while the 
MC renormalized masses and 4-point couplings are mfn and Xfn, respectively, as defined in the 
text and their RG counterparts are m and A. The MC and L-RG calculations use the same value 
of Nlat. 
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